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ABSTRACT 


The Plasma Science Experiment (PLS) on the Voyager spacecraft have provided a wealth of 
data on tlie plasma ions and electrons in tlie interplanetary medium and the magnetospheres of the 
giant planets Jupiter and Saturn. This report presents a description of the analysis used to obtain 
electron parameters (density, temperature, etc.) from the PLS electron measurements which cover 
the energy range iToni 10 eV to 5950 eV, The electron sensor (D cup) and its transmission cliarac- 
teristics are described. A derivation of the fundamental analytical expression of the reduced distri- 
bution function is given. This is followed by discussion siiowing how the electron distribution 
function fe, used in the moment integrations, can be derived from Fe- Positive ions produce a cor- 
rection current (Ion feedtliroiigh) to the measured electron current, wliich can be important to the 
measurements of the suprathernial electron component. In the case of Saturn, we show that this 
correction current, which can either add to or subtract from the measured electron current, is less 
than 207r of the measured signal at all times. Though not shown here, these feedthrough correc- 
tions are very important during the Voyager 1 Inbound pass through lo’s plasma torus. We then 
briefly comment about the corrections introduced by spacecraft charging to the Saturn encounter 
data, which can be important in regions of high density and shadow when the spacecraft can be- 
come negatively charged. 
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PLASMA ELECTRON ANALYSIS; VOYAGER PLASMA SCIENCE EXPERIMENT 

by E. C. sutler, Jr. 

SECTION 1 
INTRODUCTION 

The Plasma Science Experiment (PLS) on Voyager is a collection of potential modulated 
Faraday cups which make both positive ion and electron measurements covering the energy per 
charge range from 10 eV to 5950 eV, The PLS instrumentation has successfully measured the plas- 
ma (ions and electrons) in the interplanetary medium (Sittler and Scudder ( 1 980), Sittler et ai. 
(1981a), Belcher et al. (1981), Gazis and Lazarus (1982), Lazarus and Gazis (1983)), and the mag- 
netospheres of Jupiter (Bridge et al. (1979a,b), Scudder et al. (1981), McNutt et al. (1981), Bagenal 
and Sullivan (1981) and Saturn (Bridge et al. (1981, 1982), Sittler etal.(1981b),Hartle et al. (1982), 
Eviator et al., (1982, 1983), Goertz (1983), Sittler et al., (1983), and Lazarus and McNutt (1983)). 
The Voyager 2 PLS instrument has also played an important role in the detection of Jupiter’s mag- 
netotail beyond the orbit of Saturn (Kurth et al. (1981, 1982), Scarf et ai. (1981, 1983), Lepping 
et al. (1982, 1983) and Desch (1983)). In this report we present a fairly detailed description of 
the electron analysis which has produced the electron parameters (density, temperature, etc.) used 
in many of the above studies. The analysis described herein is most descriptive of that used for the 
most recent publication by Sittler et al. (1983) on plasma electrons in Saturn’s magnetosphere and 
is somewhat different from that used and briefly described in Scudder et al. (1981) for Jupiter. 

The ion analysis, which is fairly straiglit forward in the solar wind (Bridge et al. (1977), Belcher 
et al. (1981)) but can be considerably more difficult in the magnetospheres of Jupiter (McNutt 
et al., (1981) and Saturn (Lazarus and McNutt (1983)) will not be discussed here. 

The paper is broken up into 8 sections with the introduction given in Section 1 and a de- 
scription of the instrument and its operation given in Section 2 (a more complete description is 
given in Bridge et al. (1977)). Section 3 gives a schematic description of the D cup and its transmis- 
sion characteristics, followed in Section 4 by a formal derivation of the analytical expression used 
in our fits to the measured electron spectra; all our parameter estimations, either directly or 
indirectly, are derived from this expression. We then mention in Section 5 the fitting procedure used 
in the analysis, which plays an important role in the moment integrations. In Section 6 we expand 
upon the discussion in Scudder et al. (1981) concerning the moment estimation of electron parame- 
ters. Section 7 describes the effect of ion feedthrough coirections upon the electron measurements, 
which in the case of Saturn, introduce a minor correction to the observed suprathermal component, 
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Finally, in Section 8 we comment about the effects of spacecraft charging upon the electron mea- 
surements. 


SECTION 2 

INSTRUMENT DESCRIPTION AND OPERATION 

The PLS instrument shown in Figure 1 is composed of four potential modulated Faraday 
cups denoted by the letters A, B, C and D. The three main sensors A, B and C make only positive 
ion measurements and, except for rare brief spacecraft maneuvers, are always pointed nearly along 
the spacecraft-Sun line. The side sensor or D cup makes both positive ion and electron measure- 
ments and is normally oriented nearly at right angles to the solar direction. The angular response of 
the side sensor is cylindrically symmetric about its look direction, and provides a field of view with 
conical half angle 'v 30° (FWHM) about its normal. As shown in Section 3, the D sensor makes 
differential contiguous measurements of the electron distribution function along tlie sensor normal. 
Because tiie electron thermal speeds are much larger than Oow speeds of the plasma, electron mea- 
surements are not very sensitive to sensor orientation, unless there are large pressure anisotropies, 
Since the instrument angular field of view is fairly broad, uncertainties due to pressure anisotropies 
arc not expected to have an important effect upon the analysis. 

For cold ions in the solar wind and magnetospheres of Jupiter and Saturn the mach numbers 
are usually greater than one and the measured currents are sensitive to sensor orientation. During 
the cruise phase of the mission only the main sensors are sensitive to the supersonic ion component 
of the solar wind, while the side sensor provides a one-dimensional view of the electrons at nearly 
right angles to tire flow direction. During the encounters with the giant planets, Jupiter and Saturn, 
the D sensor was aligned to respond to the azimuthally flowing cold ions as much as possible. For 
most of the inbound portion of the encounter trajectories the D sensor was aligned to view the cold 
ions; while, except for the Voyager 1 Saturn encounter, the D sensor was not favorably aligned dur- 
ing the outbound passes. During the planetary encounters, the D sensor alignment was such that 
the center of its field of view generally looked at electrons with pitch angles between 45° and 135°. 

Except for the Cruise 1 phase of the mission when electron measurements are made every 
1 2 seconds and only the El mode is sampled, the side sensor completes a measurement cycle in 96 
seconds during which it passes through the mode sequence M, El , L and E2. M and L are the high 
and low resolution positive ion modes, respectively, while El and E2 (see Figure 4) are the low and 
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Jiigli energy electron motles, respectively. The energy range for HI is 10 eV to 140 eV while lor H2 
it is 140 eV to 5950 eV. Hach electron mode is composed of 16 contiguously spaced energy chan- 
nels; for H2 mode only the upper 1 2 channels are used (lower 4 channels not useful since suppres- 
sor is biased at 95 volts; see discussion in Section 3). The channels for HI are nearly equally spaced 
in energy (.099 ( *iK/H ( .37), while for Ii2 they are more logaritiimically spaced (-1H/H = .29). The 
sampling time for both energy modes is 3,84 seconds (0.96 second for Cruise 1 ) anci HI and E2 
modes are separated in time by 45 seconds. This large time gap between low and high electron 
ejiergy measurements can result in discontinuous changes in the composite energy spectrum across 
the 140 eV boundary joining the two energy modes. Fortunately, this happens only rarely, and the 
cold and hot components characterizing the electron distribution function within Juplte*'*® or Saturn’s 
magnetosphere are usually measured by the low and higli energy modes, respectively. The ion and 
electron measurements are never made simultaneously (the shortest time difference between ion 
and electron spectra is 25 seconds), which may lead to time aliasing problems whenever intercom- 
parisons between ion and electron measLirements are made. 


SECTION 3 

D SENSOR DESCRIPTION AND TRANSMISSION CHARACTERISTICS 


The D cup oi- side sensor is schematically displayed in Figure 2, It has a cylindrical geometry 
with ejitrance aperture at one end anci collector plate at the other end, and numerous grid meshes 
in between, The orientation of the side sensor normal relative to the spacecraft coordinate system 
is shown at the bottom of Figure 2. VVe have defined the sensor normal such that it points into 
the sensor and is thus oriented opposite to tiie sensor look direction. For added information about 
potential modulated Faraday cups we suggest reading tiie review article by Vasyliunas (1971) which 
gives an in-depth discussion of the use of potential modulated Faraday cups for space applications. 


The Faraday cup sets up a one-dimensional potential barrier, aligned along the sensor nor- 
mal, between the modulator grid (gricl3) and tiie ground grid (grid 2) shown in Figure 2. This bar- 
rier is only effective for those charged particles having a charge of the same sign as the modulator 
voltage Vjyj relative to ground potential. The dc voltage Vj^ defines the energy or speed channel at 
wiiich electrons are sampled; the superimposed 400 Hz square wave voltage, shown schematically in 
Figure 2 with amplitude ^ Vjyj, defines the energy or speed channel window size, Note that Vjyj is 
negative for electron measurements. The instantaneous current received by the collector is the inte- 
grated flux of electrons with velocity component aligned along sensor normal such that 
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with 


( 1 ) 


Vm = Vm ± A Vm/2 

and q i= ~e for electrons, A current with both a dc and 400 Hz component flows into the collector, 
Since the amplifiers are ac coupled to the collector, only those electrons satisfying the condition 


»'n ^ *^nj+ 


(3) 


where 


= Q(^Mj ± ^Vmj/2) (4) 

for speed channel j are sampled, This signal is then amplified, phase detected, and integrated before 
transmission. The role of the intermediary grids (4, 5 and 6) is mainly to reduce the capacitive coup- 
ling from modulator grid to collector plate. The suppressor grids (7 or 8 depending on instrument 
mode) main purpose is to return secondary electrons emitted by the collector back to the collec- 
tor. For electron measurements in the El and E2 modes, the suppressor voltage Vg is --8 volts 
and ~95 volts, respectively. In the normal suppressor grid configuration (grid 8 is grounded and 
grid 7 is biased at voltage Vg) the suppressor is not very effective in returning electrons back to the 
collector. Electron measurements are predominantly made in the normal grid configuration be- 
cause the instrument is considerably quieter under these conditions (Lazarus, private communica- 
tion), Preliminary estimates indicate that secondary electron corrections are only 1 0% for the 
thermal electrons (secondary electron yields are low), while *>^30% corrections are expected for 
suprathermals. At present these corrections have not yet been incorporated into the electron analy- 
sis. Because electrons with energies less than (Vg) cannot penetrate the potential barrier set up the 
suppressor grid, measurements are confined to energies greater than (Vg). It is for this reason the 
lower four E2 channels are not useable. We note that the D cup must be oriented at more than 75® 
from the solar direction, otherwise UV light striking the modulator grid will cause photoemission 
from the modulator grid, producing a contaminating signal that can swamp the measured current 
due to the plasma electrons. This problem was only intermittently present during the early phases 
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of the cruise mission before the Jupiter encounters. Because the spacecraft is stabilized w.r.t, 3 
axes, the electron measurements yield only a one-dimensional view of the electron distribution 
function, 

Transmission Function T(i^ ! 

In order to make a quantitative estimate of electron parameters such as the electron tem- 
perature from the measured currents, one must have an accurate determination of the phase space 
sampled by the sensor. This is given by the transmisison function T(p^ for which the index j 
specifies the speed channel. The transmissior,' function is defined to be that fraction of a monoener- 
getlc unidirectional beam of charged particles uniformly illuminating the entrance aperture which 
reach the collector plate. Because of cylindrical symmetry of the D cup the angular dependence of 
T is given solely by the angle of incidence relative to sensor normal o (see Figure 2). We introduce 
the normalized response R(l^ where 

(5) 

The constant Tj^j = 0.56 is the normal transmission of the sensor (i.e., 6 = 0°) and it is equal 
to the product of the transparencies of all grid meshes times the ratio of the shielding ring (shown in 
Figure 2) and aperture cross-sectional areas. In Figure 3 we have plotted a family of curves for the 
angular response RC’p; v^j-) as a function of 0 = Each curve corresponds to a differ- 

ent speed along sensor normal within the speed window of the lowest El speed channel (i.e., 
j;ni_< Vn < Vnl+). The angular response of R with half width 'v30°is principally caused by the 
common overlapping areas of the aperture and shielding ring projected upon the collector at angle 6 . 

The dependence upon or electron energy, which can be seen to be small, results from the refraction 
of electron trajectories as they pass through regions of nonzero electric field (e.g., between modu- 
lator and ground grids 2 and 4). Furthermore, the angular and energy dependence is essentially 
independent of speed channel. 

Electrons typically have large thermal speeds cog > 1000 km/s relative to plasma flow speed 
V < 600 km/s. It follows that the angular width of the electron distribution function fg seen by the 
sensor at the one thermal speed level (transverse direction, “ Wg) will be typically greater than 
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25° ill the lower speed channels. Therefore, the angular width of and R are comparable, resulting 
in a folding or convolution of instrument response with fg in the transverse dlrectlr.j relative to sen- 
sornormal. This effect is taken into account in the analysis. 


in contrast to the broadness of R relative to fg in the transverse direction, the instrument 
samples differential slices of fg along the sensor normal. To show this we present in Figure 4 a 
plot of the observed reduced distribution function Fg (similar in shape to fg) measured within 
Saturn’s extended plasma sheet at 1981 238 09:04 SCUT by the Voyager 2 spacecraft. The El and 
E2 energy or speed ranges are denoted in the figure; the histogram format is used to indicate the 
width of each speed channel with mean speed . Tlie abrupt change in channel width at 
140 eV occurs at the boundary separating the El and E2 energy ranges; the fractional window sizes 
(Av/v) for E2 are about a factor of three larger tlian that for El . The two component structure of 
the electron distribution function is clearly demonstrated by this figure. But most importantly, 
with regard to the analysis, the figure shows the differential character of the measurements in veloc- 
ity component Mathematically, this condition of differentiality, which is derived in the next 
section, is given by 




6 \C0( 


<< 1 


( 6 ) 


where 


“ *^nj+ “• *^nj- 

is the speed width of speed cliannel j, For the speed channels at whicli the thermal electron mea- 
surements are confined (El mode), Av,^j 'vSOO km/s. Referring to Eq, (6), aiong with the fact 
that cog > 1 000 km/s, one finds ej to be < 2% for all j ; therefore, the El measurements are differ- 
ential alongli'. Thus, because the speed channels are so narrow, one can in principle measure elec- 
tron temperatures less than 10^° K or 1 ev. In reality this is not always possible, because of the 
contaminating signal introduced by the suprathermal electrons near the breakpoint energy Ebi (see 
Figure 4). Furthermore, because the measurements are confined above 10 eV, tlie signal in the low- 
est El channels for temperature Tg < 1 eV will be more than 10 below the peak flux level which 
occurs at E < 1 eV. Therfore, electron densities must be sufficiently high iig > 20 to 300/cm^ 

(i.e., exact value depends on or Ijvj which are variable) with 0SC ” ^ voRs, if the cold electrons 
are to produce a signal greater than instrument threshold ~ 10 to 10^ femptoamp or instru- 
ment noise 1^ >75 femptoamps in the first few El channels. In the interplanetary medium where 
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0SC + 10 volts, a 1 eV low density (i\. < 1/cm ) electron component is detectable when" 

ever spacecraft potentials arc this higl). Tlte electron measurements in the higher energy mode K2 
can also be shown to be differential. At these liigher energies (E > 140 eV) the thermal energy of 
the electrons usually scales with electron energy H (i.e., power law in electron energy or speed); 
therefore, at I keV the thermal speeds are vl 8,000 km/s. At 1 keV the channel widths 
'V3.000 km/s, so that <?j estimated from Hq. (6) is less than ; hence, the li2 measurements are also 
differential in velocity component 1 ^,,. Witli tliese facts in mind, we will now proceed to derive the 
analytical expression for tlic reduced distribution function used in our fits. 

SECTION 4 

DERIVATION OF ANALYTICAL EXPRESSION FOR REDUCED 
DISTRIBUTION FUNCTION F^ 


To begin it is necessary to write down the general expression relating tlie measured current 
(ac component) and tlie electron distriubtion function fg(v) we are trying to determine, Keep- 
ing in mind tlie definition of the transmission function T(rs r'nj±) the fact that one can imagine 
f(,(v) to be a weighted distribution of monoenergetic unidirectional beams of particles incident upon 
tile sensor, it follows that the sampled current is given by 

. OQ OO OO 

Alj = qATN I J i\Px f f i’e(*^ )R( 

.,03 ~oo ... 


OO OO OO \ 

f / j tV,('?)R( >'nj+)''nd''n ' 


( 8 ) 


hij+ 


where A = 100 cm^ is the cross-sectional area of the entrance aperture, and ^y“- it 

can be shown (see Sittler, 1978) that these two integrals can be combined to an accuracy better 
than 1% yielding 


Alj = qATn 


OO OO ^’njf 

J dvx J di/y J 

_oo 


(9) 
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Since Mie measurements arc differential along t'n it Is useful to use the mean value theorem 
for the Vn integration, Doing this Eq. (9) becomes 


00 00 

Alj « qAT^ p nj / J (10) 

_00 — .(XJ 


where we have introduced the mean normalized response 


vnj+ 

1 "ni J R( ^1 >'ni-) •'nO-n 

•^nj- 

The speed has some value residing between Vnj± where 


( 11 ) 





1 

2 


(^’nj+ + >^nj-) 


( 12 ) 


is the mean speed for the jth speed channel. If the speed windows are sufficiently narrow, then one 
can further simplify Eq. (10) by setting Vnj* * order to estimate the error in making this 

approximation we made a Taylor series expansion of fg about and substituted it into Eq. (9). 
Doing this, assuming a step function for R in (good approximation, see Figure 3), using the 
mean response R in place of R, and using a convected-Maxweliian for fg one gets 


Alj = qATNf' nj A 


oo 00 

J ‘‘"x / dvyfeCvx, »^y, ‘^nj)R(»'l^/%j^5‘'nj J 

„oo —00 


(13) 



where V is the mean vector velocity of the electrons relative to the sensor coordinate system. Since 
Vn < 600 km/s, and > 2000 km/s, the term ^nl^nj i dropped, The result- 

ing correction term is Identical to that in Eq. (6) where it was shown in general to be less than 2%. 
The smallness of this correction term means the measurements are differential in p^ and the ap- 
proximation setting VjTj* ” has an accuracy better than 2%. Then by noting the fact that R is 
essentially identical for all speed channels, we get the following general expression for the measured 
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M 00 

“ n^'^N*'nJ^*'nj J* J ^eC»'x,»'y/jij) R(»^j2/»^nj^) (14) 

* *■■0^ mmOO 

with an accuracy better than 2%. 

In order to write down an analytical expression for Eq. (14), when for instance fg is a bi • 
Maxwellian, we fit the following sum of Gaussians to R, 

Ai e““^ (15) 

The result of this fit is shown in Figure 5, where A\ * 0.929, aj2 ■ l.l i, A 2 * 0,0813, and 
a 2 ^ = 19.09, As can be seen this function gives a v?ry accurate description of R , Substituting 
Eq. (IS) into (14) and integrating one gets 


Mj = 


qATN^nj^Aanj 


ri ■ 

® Gjk 


k=l 


n. 


3/2 2 ® 


-Wi^nj -V„)2(l-eBb„2) 


where Gjk is the integrated response 


(16) 


°jk = 


Ak e 


3k' 


V 




a 2tj 2 

^yj 


l^fni ) + ^^n r 0 + eflby^) + 


V(^j2/3i(l +eBbx^) + ak2) \ + ^gby^) + a^^) 

and fe, assumed to be a bi-Maxwellian, is given by the following expression 


(17) 


fe(^) = 


,3/2 „ 


■0Ji(^l - ^)2 + A/ir (V||- V,| )2 ) 


II ^1' 


(18) 
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The parameters' used in Kq. (16) are defined as follows! 0 ) /J|i,x ® l/w|f 
(2) 1/2 nig W|| where T|| j_ is the electron temperature parallel and perpendicular to B, 

respectively, (3) * iip/iJjL is the anisotropy parameter where ^/3 « pu - , (4) Ujj^ * (Vx/e ,^jj2^ 

Uxj'* * tVx/ e „j)- and Uyj2 » (Vy/ Vjjj)- (5) b “ is a unit vector parallel to , and ig 
is the electron density. If we assume isotropy (good approximation for thermal electrons) 
Kqs. (10) and (17) reduce to 


Alj * [qATjsjn’e,^j^Ae„J 


k«l 


'(aj{“Wx^)/(a|^2 + j/j,j2/wg2y 
“■-.V—-- ^ 

Uk*- + Vnj‘'/wg-5 




Vn)‘ 




nj vn 

?r— 

Wa** 


(19) 


where Wg is the electron thermal speed and “ Vx“/wg^. For purposes of simplicity we have 

approximated the two Gaussian fit Hq. (!5) with a single Gaussian R ~ e 
the dots in Figure 5 with a- = 1.35. By doing this liq. (19) reduces to 


A lj = [qATN7r7„j3Av,ij ] 


(a-Wx“)/(a2 +P',^j2/wg2) 


a- + F„j-/wg- 


(?nj ™ V,,)- 




3 


(20) 


which lias an accuracy better than 5‘,i for all speed channels. Finally, dividing both sides of Eq, (20) 
by qA Tjq Af^,^j one gets the reduced distribution function 


lie 


j^3/2y 3 


(?nj ~ Vn)“/(wg“) 


(a2vVx3)/{a2 + r,^j2/^.^,2) 


T/ T. 


a2+Tnj-/wg2 


( 21 ) 


This expression can be generalized to include the energy shift correction introduced by the 
spacecraft potential 0g(^ where most of the energy change is assumed to occur along the direction 
(sec Sittlcr, 1978). Doing this, Eq, (21) has the more general form 


F..i = 


Jig “ i‘SC" “ Vji)"/ ^e^ 


;r3/2u,e3 


(a2Wx2)/(a2+F„j2/wc2-j 


a2 +?nj3/wg2 


( 22 ) 


whej-c 1/2 nig usc“ - e0sC‘'ii4 Vj^j 2 > ugc 
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The reduction of the integral expression in liq. (8) to the analytical expression in Eq. (22) 
has achieved a treinendoiis simplification of tlie analysis. Instead of liaving to evaluate two 3- 
dimensional integrations for each measurement point in the spectrum during, for example, a model 
fit to thespectri'm, one has only to compute an analytical expression for each measurement point. 
Look-up ta^ie techniques could be used to enliance the computation speed of the analysis based on 
Eq, (8), but this method, whicli tends to be less llexible and not as intuitive, would still be slower 
and require more memoiy in a computer. As shown in the next section, Eq. (22) also allows one to 
directly extract from tlie measurements. In the case where the electron distribution function is 
well reprc.sented by a convected-Maxwellian, estimates of the vlectron density n^ and temperature 
can be obtained by fitting the riglit-liand side of Eq. (22) to the observed reduced distribution 
function 


A I 


F . = 


(23) 


In reality, one cannot obtain an estimate of n^, witliout having some means for estimating the space- 
craft potejitial In tlie case of tiie cruise analysis, the ion analysis was used to give an estimate 
of the electron density and charge neutrality was enforced. This approach, which gave a direct 
measure of the spacecraft potential, was used to construct tire return current relation presented in 
Scudder et al. (1981). This return current relation was then used as a constraint equation in tlie 
analysis which gave a self-consistent estimate of both the spacecraft potential and electron density. 
The use of Eq, (22) is also limited by the fact that is not described by a single Maxwellian for 
the full energy range of the plasma instrument. As discussed extensively in Scudder et al (1981) 
and Sittlcrctal, (1983), the distribution function is characterized by a cold and hot component, 
wliere the cold component is Maxwellian in form, while the hot component is non-Maxwellian in 
its energy dependence (see Figure 4). To attack tins problem we have fit a sum of three Maxwel- 
lians to the reduced distribution function F^, where Eq. (22) is used for each Maxwellian compon- 
ent. For the spectrum in Figure 4 we have superimposed the result of this three Maxwellian fit to 
Fg using the x symbol to denote the fit value at each speed channel. A single Maxwellian describes 
the cold component quite well; two Maxwellians give a fair description of the suprathermals, for 
which a power law would probably give a better fit. As discussed in the next section, the Maxwel- 
lian fits to the hot component are only performed for thy purpose of estimating the distribution 
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function fg from Fgi fg is then used in the moment integrations to compute moment estimates of 
ng and Tg. 

The details of the analysis procedures used to perform these fits can become quite complex 
because of the wide range of plasma regimes sampled; for example, one must guard against the 
possible contamination by photoelectrons from the photoelectron sheath surrounding the space- 
craft, In most simple terms, the analysis is performed by first fitting a Maxwellian to the cold elec- 
trons in the first few speed channels. In the magnetosheath and magnetosphere regions, the flow 
velocity has been set equal to zero, and in the solar wind the bulk speed from the ion analysis is 
used. Once this has been done, the breakpoint energy Egj is computed, which is defined to be 
the energy at which Fg is twice the value of the Maxwellian fit ^MC to the cold component F^ 

(see Figure 4). We then fit a Maxwellian to *^H1 “ *^e ” ^MC in the third to seventh channels above 
Eg| , The Maxwellian Fj^^jj fit to Fj^j is then subtracted from Fg in the energy channel below 
Egj (i.e„ F^ ■ Fg - Fj^i^j). We then fit a Maxwellian to this revised F^ and then repeated the 
same above procedure used to compute Egj and obtain the Maxwellian fit to Fp|j. A similar pro- 
cedure is used to get a Maxwellian representation of the hot component Fj ^2 above the second break- 
point energy Eg 2 > As can be seen from Figure 4 this procedure yields a good description of the 
observed reduced distribution function Fg. The Maxwellian fits to the cold electron component arf’ 
usually very good, on average the goodness of fit parameter Xg^ “ 5, The for the Maxwellian fits 
to the the suprathermals are larger, being typically greater than 10, but usually give a good estimate 
of the thermal width of the spectrum over the energy range each Maxwellian component is fit to. 


SECTION 6 
MOMENTS 

As discussed in Scudder et al. (1981) it is desirabie to have a quantitative estimate of the 
total density and mean random energy of the electron gas as a whole, since it is these quantities 
which appear in the fluid equations. The fitting of a Maxwellian to the distribution function is 
only useful when it gives a good description of the particle population. Clearly, a Maxwellian 
does not give a good representation of fg for the full energy range of the PLS instrument, but does 
give a, good description of the sub-interval of energies occupied by the thermal electrons. The sup- 
rathermal electrons are definitely non-Maxwellian; as shown in Section 5, more than one 
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Maxwellian is needed to describe them. A power law or kappa distribution could be used and would 
probably give a good but model dependent description of their energy dependence, The reduced 
distribution function Fg, contrary to f^ which has a direct connection with theory, contains the 
effects of the instrument function. Furthermore, once fg is known, it is a relatively simple exercise 
to integrate over fg and obtain model independent estimates of the density and temperature of the 
gas as a whole. 

As noted in Scudder et al, (1981) the distribution function fg can be retrieved from the 
measurements, so Jong as the speed window me is small compared to tlie local thermal speed 
of fg. By local thermal speed we mean the magnitude of the square root of the local logarithmic 
derivative of fg with respect to electron speed. Careful inspection of Eq, (22) reveals tliat one can 
estimate fg from Fg if the distribution function is Maxwellian and tiie thermal speed o>g is known. 
Therefore 


fgj “ tg(o,o, i^y) - I-gj(a- + V|^j‘'/cOg“) ^24] 

where we have dropped die exponential term with argument proportional to W^-, since for the 
electrons < < 1 , One can carry this approacli one step further by saying that f^j can be esti- 
mated with Eq, (24) even if fg is not a Maxwellian at by using in place of Wg in Eq, (24) the 
locally determined thermal speed wgj = o>g (F,ij) of fg at speed F]^, Note that the correction factor 
a~ + (VnjVcjg)^ is an indication of the angular width of fg relative to the sensor response R at 
speed F,ij, Therefore, by estimating the thermal width of fg along and assuming this 
gives a good representation of the thermal width of fg in the tramiverse direction (isotropy assum- 
ption), Eq. (24) should give an accurate estimate of fg atF^j. In order to make use of Eq, (24) we 
must have an estimate of cogj, Two approaches have been used. In the first approach we make a 
two point estimate of cogj from Fgj and Fgj+i and using Eq, (22), In the second approach we make 
use of the three Maxwellian fits to Fg to estimate Wgj. Centered on the breakpoint energies Eg] 
and Eg2 where the thermal speed Wg has a strong speed dependence, we have linearly interpolated 
across the breakpoint energies Egi and Eg2 to obtain an estimate of cogj in these regions. Both 
methods yield similar results whenever the signal is well above instrument noise or temporal fluctu- 
ations are absent. In the lower density regions or whenever rapid temporal variations are present 
in the spectra, the first method has problems because the noise or temporal fluctuations introduce 
localized changes in the slope of Fg which then introduce large errors in the estimation of ajgj. 
Therefore, we have used the second approach which is less sensitive to noise or temporal fluctua- 
tions since the thermal speed is estimated over a broader e. y range and more points are used in 


13 


ORIGINAL PAGE IS 
OF POOR QUALITY 


estimating the local slope of f^ (i.e„ thermal speed estimation less sensitive to single point fluctu- 
ations in the energy spectra). One means of checking this procedure for estimating fg, is to see how 
well the Maxwellian fits in Fg vs “phase space” compare in fg vs Pl’^^se space. If the conver- 
sion ha« been done properly, and the fit to Fg is of a good quality, then the same fit parameters 
should yield a good match to fg where 

f j s e ~ (v^nj“ " usc“ V„)-/wg2 

rr*’/- jg-^ 

is used for each Maxwellian component. Figure 6, wJiich is the same spectrum plotted in Figure 
4 except that fg is plotted instead of Fg, shows the Maxwellian fits give a good match to botli Fg 
and fg. For all spectra this comparison has been made we found a similarly good match. 


Once the distribution function fg is known, one can proceed to nmke moment estimates of 
the total density iig and temperature Tg as outlined in Scudder et al. (1981 ), For reference pur- 
poses we write tlie integral expressions for iig and Tg, respectively, 


llg S 4 7T 


OO 

J ‘V^d)''d“ tl^d 

0 


(26) 


CX) 

^ J ^V'^d) 


1 ^ 


»^d/(-f'^e 1<B) 


(27) 


where 



2e0sc/”ie 


Yz 


(28) 


Is the electron speed outside tlie pliotoelectron sheath surrounding the spacecraft. The speed uq 
is the observed electron speed at the spacecraft surface, e is the unit electric charge, nig is the elec- 
tron mass, kg is Boltzmann’s constant, and <P^q is tlie spacecraft potential. The above integrals 
assume isotropy in fg where we have taken advantage of the subsonic character of electrons and 
assumed corrections introduced by pressure anisotropies to be small. As a first approximation, 
we have only considered the energy shift correction caused by the spacecraft potential in Eqs. (26) 
to (28). In Sittler et al. ( 1 983) the upper limit used in the integrals 26 and 27 is set by the maxi- 
mum channel electron fluxes arc observed above Ij[.j or 4’*Tjsj (Ijsf is the theoretically estimated 
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lower limit of the instrument noise; in reality the instrument noise can be greater tlian ij^), sucli 
that tlie maximum energy of tiic integral can be no more than 6 kcV. Because the bulk of the dis- 
tribution is below 6 keV witliin Saturn’s magnetosphere» tlie upper limit for the integrations is 
essentially infinity. In Scudderetal. (1981) we extrapolated the integrals to infinity, since the 
suprathcrmal electrons in Jupiter’s magnetosphere were hotter ~ few keV) and noise in the 
higher channels less prevalent. 

We now amplify the point made in Scudder et al. (1981) concerning the extrapolation 
down to zero energy. Because the measurements are confined above 10 eV, one must use the esti- 
mated temperature for the cold electrons to extrapolate f^ down to zero energy in order to com- 
plete the moment integrations. As argued in Scudder et al., this portion of tlie integral will not 
make a significant contribution to the integral if f^ does not vary too much below 10 eV or deviate 
significantly from that indicated above 10 eV. The first condition is violated whenever < 10 eV 
(density dominated by cold electrons) since the integrals will peak near or below 10 cV. Further- 
more, as the electron temperatures become lower (T^, < 3 eV) most cold electrons have energies 
< 10 eV, and fluxes above 10 eV can be dominated by supratlicrmals so that it becomes more and 
more difficult to make an accurate estimate of If supratliermals are absent, the cold electron 
fluxes can be below instm.ment threshold and/or instrument noise. At those low temperatures, 
the error in estimating may be no more tlian a factor of two, but since the spectrum is so 
steep below 10 cV, large erorrs in the density can result. This accounts for some of the data gaps 
in the density plots in Sittler et al. (1983), while no data gap is present in the temperature plots; 
in the same plots other regions have been omitted altogether. Eventually, by using the ion and elec- 
tron data sets together, the successful analysis in these cooler regions in Saturn’s magnetosphere 
is anticipated. 

Another complication arises whenever the temperature of the thermal electrons Tj. ex- 
ceeds '\^20 eV, such that a turnover in the spectrum in the lower energy channels tends to occur. 
That is, the predicted phase density for energies below 20 eV from a Maxwellian fit to energy 
channels greater than 20 eV, exceeds that observed below 20 eV. This low energy turnover in the 
spectrum has two possible explanations: (1) we arc seeing an effect similar to that seen in the 
Earth’s magnetosheath where the eiectron distribution function is observed to have a flat 
topped shape at lower energies; ion acoustic waves have been proposed as the means for providing 
this distortion of fg (Dum et ah, 1974); or (2) it is an instrumental effect caused by the emission 
of secondary electrons from the modulator grid. This latter mechanism, first proposed by Vasyii- 
unas (1971) is caused by the change in energy of primary electrons incident upon the modulator 
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grid as the modulation voltage Is varied; by changing their energy the secondary electron yield, 
which is a function of the primary electron energy, will also change. Then if the flux of electrons 
is greater at higlier speed channels, a significant correction to the lower electron channel is feas- 
ible, where this correction current will tend to reduce the measured signal in the lower electron 
channels (i.e,, below 100 eV the secondary yield .is an increasing function of primary electron 
energy). Since the flux maximizes at the electron thermal speed, this effect, if important, will 
occur whenever Tc ^20 eV and will affect only those speed channels below cop. 

Detailed calculations of this effect are being planned for the future. For now we have cir- 
cumvented this problem by confining our analysis to speed channels above this turnover in the spec- 
trum. If the first explanation applies, then as noted above, we will make only a modest overesti- 
mate of the electron density by extrapolating to zero energy the Maxwellian fit to speed channels 
above the low energy turnover in the spectrum. If it is an instrumental effect, then we make no 
error by extrapolating to zero energy the Maxwellian fit to speed channels above this turnover in 
the spectrum. 

Since the suprathermal electrons in Saturn’s magnetosphere are non-Maxwellian, we have 
used moment estimates of their density and temperature in our discussions about them in Sittler 
et al. (1983). As in Scudder et al, (1981) the lower limit of the integration is set by the 
breakpoint energy Eg j; the upper limit is the same as that used for the density and temperature 
determinations in Sittler et al. Note that this method will tend to underestimate the density iij^jand 
pressure Pj^ of the hot electrons since we arc ignoring their phase density below Eg| . This was done to 
avoid numerical problems which occur if we equate the fractional density to the difference 
n„ - n„ where n„is the fit density of the cold electrons. These numerical problems will occur when- 
ever nj^/Uj, S 0.05 or when the fit to the cold electrons is poor (which will happen when El mode 
fluxes arc near instrument noise). This method is also desirable because it allows one to study more 
clearly the variations of the more energetic suprathermal electrons above the breakpoint energy Eg j . 
Before the integrations are performed we first subtract the cold component from f^ by using the 
Maxwellian fit to the cold electrons. Finally, whenever the measured current is below ItH 
the distribution function f. is set to equal zero. This approach, which is also used for the n and T„ 

v< C w 

integrations, eliminates the possibility of over estimating the integrals because of noise, while tending 
to under estimate the integrals because signal less than lyp£ or may be present. 
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SECTION 7 

ION FEEDTHROUGH CORRECTIONS 

In Scudder et al. (1981) we made a point that because the instrument records and teleme- 
ters only an ac synchronously detected signal, dc background corrections are not a problem. For 
the high energy particle population this is a correct statement, But low energy positive ions, which 
have energies per charges < 6 keV comparable in magnitude to the potential applied to the modula- 
tor grid of the D cup, will produce an in phase ac signal when the instrument is in the electron 
mode that may or may not introduce an important correction to the measured electron currents. 
This feedthrough current originally pointed out by Vasyliunas (1971) as a possible correction 
current, is caused by slight changes in the ion tra>ctories (see Figure 2) as they pass through 
regions of nonzero electric field within the sensor (e.g., the modulator grid), When the voltage on 
the modulator grid is varied at 400 Hz, ions with same energy, angle of incidence, and point of entry 
will experience slightly different deflections in their trajectories. This will cause some of them to 
either hit or miss the collector plate (which is ac coupled to the amplifier network) at a 400 Hz rate. 
The end result, as noted above, is an in phase ac signal to the amplifier network, The sign of this ion 
feedthrough current is such that it will add to the electron current Ig (e.g., lQj 3 g ~ ie ^fed 
where lQg§ is the observed current). 

Using the simulation program developed by Vasyliunas (see Hartle et al., 1 982, in which elec- 
tron feedthrough corrections to the ion measurements were computed) and predicted ion parameters 
at the spacecraft position from scale height model calculations, we computed estimates of ion feed- 
through currents for all electron energy channels througliout a major part of the Saturn encounter 
period for both Voyager encounters. The results displayed in Figures 7 and 8 show that the ob- 
served electron currents due to the thermal electrons are more than two orders of magnitude 
greater than these correction currents, an expected result. For the suprathermal electrons, for 
which these corrections have the potential of being important, the feedthrough currents are found 
to be no more than 20% of the observed signal. As demonstrated in Sittler et al. (1983), the 
suprathermals contribute about 10% to 20% of the total electron density. It then follows that 
these feedthrougli currents are expec :ed to be important whenever Uj.j/ng ^ 1%, a situation which 
is known to occur during the Voyager 1 inbound pass through lo’s plasma torus at Jupiter. Since 
these corrections are only minor at Saturn, we have not included them in the analysis presented in 
Sittler etal, Whenever the ion analysis is completed, these correction currents will be incorporated 
into the final electron analysis that is to be part of a future ion and electron paper for Saturn. 
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You will note tlvat there are occasions at the higlier energy channels (channels 25 to 28) 
when the feedthrough current exceeds the measured current Iqqs* At first thought this is 
difficult to understand since Iqbs ® le ^tcd should always be greater than Ife^i. Focusing our 
attention on the Voyager 1 data from 0100 to 0400 on day 318 (similar arguments can be made 
for the Voyager 2 data near the end of day 237), we see that the observed spectra display a strongly 
energy dependent behavior at higher energies and that signal is confined near or below signal thresh- 
old at the higher energy channels (bite-out signatures discussed in Sittler et al, (1981b, 1983)), This 
should be contrasted with the relatively flat spectrum characteristic of the feedthrough current 
*fod Pi'otlticed by the cold ions. VVe further note that the density profile of the cold ions produc- 
ing displays a gradual rise during this same time period {'Lazarus and McNutt, 1983), while 
Iqbs electrons can be quite variable. So the variations in the electron current, which 

can be quite variable and display a significantly stronger energy dependence than If^j, cannot be 
caused by density variations in the ions, but rather are caused by real changes in the hot electron 
fluxes. The above discrepancy can be understood in terms of a feedthrough current with 
sign opposite to that of , that is produced by ion impact upon the modulator grid with subse- 
quent secondary electron emission from it. This feedthrough current I(^^i , based upon preliminary 
estimates discussed below, is generally small compared to except at the high energy channels. 
Therefore, at these higher channels the feedthrough current is expected to be small and possibly 
negative. 


The mechanism producing this negative feedthrough current is similar to that discussed in 
Section 6, except the incident particles are ions and not electrons. Because the positive ions have 
polarities opposite to that applied to the modulator grid, all the incident ions are accelerated to- 
ward the modulator grid between grids 2 and 3 (see Figure 2). By varying the voltage applied to 
the modulator grid, the energy at which the ions strike it is varied. Then if the secondary yield 
is energy dependent, the secondary electron current is also varied, This feedthrough current will 
be out of phase (negative) with Ig if the secondary yield Yj is an increasing function of ion energy 
Ej, For Ej < 1 keVthe secondary yield Y| is <0.1 and essentially independent of Ej (Whipple, 
1965). Therefore, for those energy channels the ions strike the modulator grid with Ej < 1 keV, the 
feedthrough current will be near zero and unimportant.. In the case of tungsten (material 
grid wires made from) and protons being the incident ion, the yield Yj is an increasing function of 
ion energy 1 keV < Ej < 100 keV (Whipple, 1965). A similar energy dependence is expected for 
0"^ ions. Therefore, the feedthrougli current will be negative and reduce the observed current 
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Jobs * Jc J|cd + Jfed * an incjx'iising function of ion energy and tiie energy win- 

dow size qj ^ Vjyjj increases with Vj^j (E2 mode), I will be most important at the higher energy 
channels for which the ions are accelerated up to energies near 6 keV before striking the modulator 
grid and secondary yields are approaching 1 (Yj = 0.5 at Ej » 5 keV for protons). 


In Figure 9, we luive plotted IqbS' Jfed’ J fed Jfed ” * fed J fed versus channel num- 
ber for plasma data measured near the Voyager 1 ring plane crossing at Saturn. The ion parameters 
used to compute and I were supplied by R. L. McNutt (private communication) for the 
ion spectrum measured nearly time concident with the electron spectrum. Except for channels 
25 to 28, IfJti !s less than 209^ of the observed current Iqbs* J fed uegligible below channel 22, 
above which it rises rapidly and exceeds Ifed at channel 28, the uppermost channel. When we combine 
^fed ed' feedthrough current is no more tiuui 109? of IqbS’ observed 

variations in Iqbs between 0100 and 0400 on day 318 for Voyager 1 and 1930 to 2400 on day 
237 for Voyager 2 arc most assuredly due to real changes in suprathermal electron fluxes. 


SECTION 8 

SPACECRAFT CHARGING EFFECTS AT SATURN 

The electron analysis for the Saturn and Jupiter data sets uses an empirical return current 
law to estimate (p^Q', this law only considers the balance between the photoelectron current emitted 
by the spacecraft and the incident plasma electron current on to tlie spacecraft. This approach is 
sufficiently accurate for electron parameter determinations, so long as the incident plasma electron 
current outside the sheatli is not high enough to enforce a negative potential on the spacecraft. This 
will happen when the plasma electron current exceeds the photoemission saturation current JsaT 
as defined in Scudder et al„ plus the incident positive ion current and secondary electron emission 
current. Tlie preliminary estimate for JgAT Scudder et al., is now estimated to be about a fac- 
tor of 3 lower, though this estimate will probably require further revision, as we did not take into 
full account the positive ion current above 6 keV and associated secondary electron emission cur- 
rent. Using this revised value for Jsat> ' ve expect the spacecraft to remain positive in Saturn’s 
magnetosphere so long as election densities remain below lO/cm^ for Tg = 10 eV. Since this 
estimate has not included the ion and secondary electron emission currents, the spacecraft could 
remain positive for even higher densities and temperatures. Except for the Voyager 1 ring plane 
crossing period at Saturn this condition is probably met. Detailed comparisons with ion data so 
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far support this conclusion, It is for this reason we feel the electron analysis of the Saturn data 
set is in good shape at this time, except for periods of shadow and high density (ne > 10/cm^) 
when the spacecraft may have become negatively charged. 
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FIGURE CAPTIONS 


Figure 1, Photograph oHlie Voyager Plasma Science lixperiment showing ihc three main scnsoi-s 
A, 13, ami C and tite side sensor 1). Contrary to tiie main sensors, tlie D cup iias tlie usual cylindri' 
cal geometry, 


Figure 2. Scliematic cross-sectional view ofD-cup. Trajectories of incident primary electron and 
emitted secondary electron are displayed. Figure shows capacitive coupling between cdiector 
plate and amplifier chain. Only current is amplified, phase detected and telemeter.'d, Voltage 
levels indicated for suppressor grid pertain to when the instrument is in El mode; when in E2 mode, 
Vg « volts. Angled and primary electron normal and transverse velocity component axes 
are shown, Tlie lower left inset shows the orientation of the sensor coordinate axis (Xs» 
relative to payload coordinate axis (X$/C' Ys/C» ^S/C^* Angles 


1‘igure 3. Normalized angular response of D-cup R(a; „) plotted versus angle of incidence d . 
Instrument is in HI mode, lowest speed channel, and modulator set at lower speed level yui... 
Each curve corresponds to electron velocity component />n(k) *i^nl- Aii>jji(K/5) for K ■ 1 to 5; 
the uppermost curve corresponds to K « 5 (t^n * «hil+^' See text for definitions. 


Figure 4. Plot of Voyager 2 electron reduced distribution function F^ plotted versus normal com- 
ponent of electron velocity ujj into D-cup, Spectrum measured in Saturn’s magnetospliere during 
outbound pass at L * 8 . Histogram format used tor data, X’s used to denote su i of tliree M.^xwel- 
lian fit to data. Cold component Fy, and hot components F^l and Fj |2 are indicated. Breakpoint 
point energies Eg] and E 132 defined in Section 5 are denoted as well as speed range covered by El 
and E2 modes, Tiie one telemetry count level (one TMO and noise level are displayed. 

Figure 5. Plot of mean normalized response plotted versus angle of incidence 

0 of primary electron. _ Mean response computed from normalized response plotted in Figure 3. 
Actual mean response R terminates at (i » ±60°, The two gaussian fit is indicated by tlie solid 
cuive which exceeds R at t? « 0° and is non-zero beyond 0 - ±60°. Circles used to indicate single 
gausian fit with a- = 1.35. 


Figure 6, Electron speed distribution function fg (circles) computed from electron reduced distri- 
bution function F^ in Figure 4, Cold component fg and iiot components fj-[i and fH2 are denoted; 
breakpoint energies Eb] and £32 are rdso Indicated, Solid curve displays sum of three Maxwellian 
fit to data in Figure 4 (see Section 5). ^ -ne telemetry count level rone TMC) and noise level are 
displayed, 


Figure 7, Plot versus time of the measured electron current lops predicted ion feedthrougli 
current for all 28 energy clnm.nels (channels 1 to 16 for El ; channels 5 to 16 for E2); currents 
are in femptoamps. The data was measured during the Voyager 1 encounter with Saturn. The 
ion parameters arc computed from a scale height model calculation; uncertainties in composition 
are not expected to produce significant errors in our estimation of If^jj • The data plotted with an 
X symbol at 0406 day 318 were computed from ion parameters determined from a detailed fit to 
the ion spectrum at this time (ion parameters provided by R. L. McNutt, private communication). 

As can be seen the model and fit ion parameters yield nearly identical results. Noise and one 
telemetry count (TMC) levels are also displayed. 


Figure 8 , Same as Figure 7 except Voyager 2 Saturn encounter data are used, 
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Figure 9. Observed electron current Ions ^solkl line) ami ion feedthrough currents IfJd (dotted 
line), Ifed (short dash) and Ifgd ■ Ifed + Ifell (long dash) plotted versus channel number for the 
Voyager I electron spectrum measured at 318 04 06 SCET; currents are in femptoamps. The feed* 
througli current ijed represents an upper estimate, but is probably correct within a factor of 2. 
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